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Abstract — In this paper, we investigate the cumulative distri- 
bution function (CDF) of the aggregate interference in carrier 
sensing multiple access/collision avoidance (CSMA/CA) networks 
measured at an arbitrary time and position. We assume that 
nodes are deployed in an infinite two-dimensional plane by 
Poisson point process (PPP) and the channel model follows the 
singular path loss function and Rayleigh fading. To find the 
effective active node density we analyze the distributed coordinate 
function (DCF) dynamics in a common sensing area and obtain 
the steady-state power distribution within a spatial disk of radius 
R/2, where R is the effective carrier sensing distance. The results 
of massive simulation using Network Simulator-2 (NS-2) show a 
high correlation with the derived CDF. 

Index Terms — Aggregate interference, CSMA/CA, DCF, Pois- 
son point process, NS-2 



I. Introduction 



A. Motivations 



Due to the inherent scarcity of frequency spectrum and 
increasing wireless traffic demands nowadays, frequency reuse 
has become an essential key technological issue associated 
with contemporary wireless communication systems. Fre- 
quency reuse intrinsically causes interference between wireless 
links in not only homogeneous but also heterogeneous systems 
using the same frequency. Accordingly, the state of the aggre- 
gate interference at an arbitrary position in the random node 
topology has become of great importance. 

In this paper, we are interested in the aggregate inter- 
ference of carrier sense multiple access/collision avoidance 
(CSMA/CA) networks. Our motivation is two-folds: The first 
motivation is associated with the prevalence of this random 
access network. WiFi that complies with IEEE 802.11 stan- 
dards offloads the traffic from conventional wireless cellular 
systems. In reality, this support from WiFi systems helps 
current cellular systems to manage the ever-increasing traffic 
without urgent and huge amounts of investments in system 
upgrades, and this is likely to last for a while. Additional 
merits of the WiFi system include the fact that it is cheap 
and easy to set up without any spectrum-licensing procedures. 
As a result of these benefits, WiFi access points (APs) and 
compliant stations such as laptops, smartphones, and even 
smart TVs are being deployed more frequently by end users 
and office administrators, as well as wireless service providers 
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in homes, offices, automobiles, and public areas. Even more 
of these devices will be deployed when advanced amendment 
standards such as IEEE 802.11 ac, ad, af, and ai are finalized 
and implemented. 

Our second motivation is the fact that it is more difficult 
to control interference in this random access network than in 
centralized systems even though the interference is crucial in 
both systems. Unlike the former system, which is designed for 
proactively controlling the interference between the commu- 
nicating entities and managing the spectrum resources in the 
central control unit, the latter system has no global control 
unit. A mechanism called intercell interference coordination 
(ICIC) [1], [2] has already been proposed to resolve the 
mutual interference in the former system and it is highly 
likely that it will be implemented in practice, although some 
amount of modification may be necessary. In contrast, WiFi 
was designed without considering interference among over- 
lapped basic service sets (OBSSs). In addition, coordination 
between multiple WiFi cells calls for policy-level negotiation 
among different WiFi service operators, which is an untenable 
issue in the technology domain. Therefore, controlling mutual 
interference in random access networks is more problematic 
than in controlled counter systems. 

For these reasons, we believe that studying the charac- 
teristics of the aggregate interference with a view toward 
using this knowledge to craft realistic solutions can solve 
the uncontrollable interference problem in the WiFi arena. 
Within this context, we investigate the aggregate interference 
in randomly deployed IEEE 802.11 distributed coordination 
function (DCF) networks (we prefer this standard name instead 
of simply WiFi so as to explicitly incorporate the PHY and 
MAC layer operations) in this paper. From an understanding 
of the distribution of the aggregate interference at the protocol 
level, we can control this interference using the relationships 
discovered among the protocol parameters. Since to date there 
has been no test at the simulator level for aggregate interfer- 
ence, we test as well as analyze the interference at the protocol 
level. It is difficult to apply the measurements from the test 
bed to statistical inference on interference power distribution 
because signals are vulnerable to many unknown factors in the 
testing circumstance and massive variations in node geometry 
are hard to be set up in the real test bed. Simulator-level 
tests can however overcome all of these weaknesses associated 
with real measurements. There is no unknown factors to affect 
aggregate power in the simulator and numerous changes in 
node geometry are possible. Consequently, our goal in this 
study is to obtain the statistical inference of the aggregate 
interference and verify the results via simulation. 



2 



Most of the work previously done in this area focused on 
ALOHA-like systems in which the aggregate interference can 
be analyzed by assuming that transmitting nodes have inde- 
pendent locations and behaviors [3], [4]. Although broadband 
cellular systems such as LTE, LTE-A, WCDMA or its femto 
cell networks can also be modeled using this transmission- 
independent behavior, this is not a realistic assumption for 
CSMA/CA networks in which a certain interference level 
always needs to be maintained in a distributed manner. In 
a network of CSMA/CA nodes, every communication entity 
first senses the ongoing transmission in the channel and then 
determines when to start transmitting. The inappropriateness 
of the independent model in such a scenario was noted in 
[5], in which the authors proposed the alternative dependent 
point process to mimic real CSMA/CA networks. However, 
this proposed point process still cannot describe the DCF op- 
eration, in which collision and idle time, as well as successful 
transmission, can occur even in the exclusion area by carrier 
sense. Compared with these previous research efforts, to the 
best of our knowledge, our work is the first to investigate 
the exact distribution of aggregate interference in CSMA/CA 
networks and to validate it by means of massive simulations. 

B. Summary of contributions and organization of this paper 
Our paper has the following notable results: 

• Owing to the possibility of concurrent transmission 
incurred by DCF operation occurring within an ex- 
clusion area, CSMA/CA random networks can only 
be modeled by the Poisson point process (PPP), not 
the dependent point process. Section [TT] addresses the 
difference between the dependent and independent point 
processes, and explains how real CSMA/CA networks 
can be dealt with. 

• We derive the effective node density reflecting 
CSMA/CA MAC layer operation. Section HITl explains 
this. This section is the core of our work since all the 
cross-layer parameters are coupled to model the network 
behavior. Further, this is used in PPP shot noise analysis, 
which is the result of the first item. 

• Aggregate interference in PPP shot noise analysis 
using our derived effective node density is verified 
using the NS-2 network simulator and Matlab sim- 
ulations. Unlike other related theoretical analyses of 
stochastic geometry researches, we verify exactly how 
much our model reflects the NS-2 simulation results. 
We also compare our results with the Matlab simulation 
on dependent point processes and thereby show that the 
PPP model with our new effective node density performs 
the best in modeling the aggregate interference. Related 
issues are outlined in Section [IV] 

Concluding remarks are given in Section[V] A compact form 
of shot noise analysis [6] follows in the Appendix. 

II. Point process for modeling random CSMA/CA 

NETWORKS 

In this section, we focus on determining which type of 
point process is suitable for modeling CSMA/CA networks. 



A generic wireless network consisting of multiple randomly 
deployed nodes can be described via a point process. In the 
point process, a mark (a scalar or a vector) can be assigned 
to each point independently, which is useful for modeling 
node-oriented properties such as transmission power, medium 
access delay. In particular, the case where the number of nodes 
in a network is Poisson-distributed and their positions at a 
given time instant are independent of each other, is adequately 
explained by means of the PPP. The method to derive the 
aggregate power emitted from points at an arbitrary position 
under independent marked PPP was previously studied as shot 
noise field, which was originally used to model the noise in 
electronic circuits in the time domain [7]. 

A. Inappropriateness of PPP for CSMA/CA modeling 

However, the PPP approach as it currently exists is insuffi- 
cient to model the CSMA/CA. The reason is that the carrier 
sensing philosophy is not reflected in it. PPP is a typical 
independent point process in which the points are deployed 
independently of each other. On the other hand, in the carrier 
sensing (CS) operation, a sensing node always senses the 
shared medium and it delays its transmission once it senses 
that the medium is busy. The result is that active nodes are 
affected by each other, which means that the process is not 
independent. 

B. Inappropriateness of the dependent point process 

Let us now consider the dependent point process as a possi- 
ble alternative. Here, the dependent point process means that 
some initially deployed points are discarded or selected by the 
metric related to the other points' marks or locations. To the 
best of our knowledge, there are two dependent point processes 
most related to the modeling of CSMA/CA networks: the 
Matern hardcore (MHC) process and the simple sequential 
inhibition (SSI)- point process. In [5], the authors are also mo- 
tivated by the inappropriateness of PPP, i.e., the independence 
of points. In the paper, they compared the aggregate power 
distributions of PPP, MHC, and SSI with simulations, and 
concluded that SSI is the most appropriate one for modeling 
CSMA networks. However, their result is not fully acceptable 
because the operation they used in the simulation was not the 
real one but a modified version of a dependent point process 
based on IEEE 802.15.4 PHY parameters. More specifically, 
they neither considered any details of practical MAC layer 
parameters nor channel characteristics. Leaving this result 
aside, their MHC and SSI are fundamentally based on the 
hard exclusion area; nevertheless, they contain ambiguities in 
the determination of this area. To further illustrate the issue, 
let us look at our simulation results. 

We simulated a realistic CSMA/CA random network using 
NS-2 in order to observe the concurrent transmission behavior. 
Fig. Q] shows our simulation topology and the effect on the 
distribution of the number of concurrent transmitting nodes. In 
the grid topology, the black dot represents the transmitter and 
the corresponding receiver is located at 5 (m) right and 5 (m) 
up away from its transmitter; the receiver is omitted from the 
figure. In the simulation, a 500 B payload was given to each 
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(a) 



Number of transmitting nodes 

(b) 



Fig. I . The number of transmitting nodes concurrently: (a) scenario and (b) 
result. 



transmitter and the traffic was saturated. The distance between 
the two nearest black dots was 50 (m), and the CS threshold 
was tuned so that the resulting CS range was 70 (m), derived 
using Equation ©. All other parameters were the same as 
those in Table [III] The large circle in the center denotes a CS 
area based on the white dot (one of the transmitters), while the 
smaller circle has a radius that is one half of the CS range. The 
relative time durations for each set of concurrent transmissions 
are depicted in the bar plot. 

In the graph, two things are of note: First, there is a time 
period in which two nodes are concurrently transmitting when 
all transmitters are even in each other's CS area (see the 4 
nodes case in the bar plot). Second, there is a time period in 
which the medium is idle in the full CS area (node 9 case). The 
first case occurs due to CS failure or collision in real situation. 
MHC and SSI are fully dependent thinning of PPP with 
the exclusion area, and they cannot model these events. The 
resulting effective node density of the dependent point process 
is likely to be lower than that of the real one. These approaches 
may work well in collision-less CSMA/CA networks where 
slot time is zero and backoff time is a continuous random 
variable (RV) [8], rather than real situations. The second case 
occurs due to the idle time from the binary exponential backoff 
(BEB) and the DCF. This waste of time resource is the intrinsic 
cost of the distributed random access MAC. In the dependent 
point process, any point having none of the other points in its 
exclusion area always survives. The resulting effective node 
density of this process is likely to be higher than that of 
the real one. As a result, the real operation of a CSMA/CA 
network has both factors having higher and lower effective 
node density than that of the dependent point process. This 
difference is from the lack of MAC layer operation modeling 
in the dependent point process. 

C. Review of PPP 

As will be shown in Section IIII-B.ll these collision events 
occur with a certain probability in real situations. This means 
that the concurrent transmission in an exclusion area oc- 
curs with some probability, not with deterministic patterns. 
This stochastic characteristic of real networks is appropriately 
modeled using the independent point process. Therefore, we 



believe that a possible way to model a CSMA/CA network is 
again to use PPP, but with a new effective active node density 
reflecting MAC layer operation. 

This is notable since recent research efforts such as [4] and 
literatures therein point out that CSMA/CA can be modeled 
with MHC point process, which is different from our conclu- 
sion. At the very least, the aggregate power at an arbitrary 
position can be elaborated more when using PPP with a new 
effective node density rather than pure PPP or MHC/SSI. 

III. Effective Active Node Density 

In the infinite plane, the nodes in a CSMCA/CA network 
are coupled with each other via the shared medium. The 
number of neighbor nodes, relative distances among nodes, 
and the current state/history of MAC and PHY/MAC/NET 
layer parameters of nodes in the local shared medium are again 
the factors affecting the states of their media. Moreover, this 
medium itself does not have a deterministic boundary but a 
random one since the channel link gain is varies continuously. 
Therefore, obtaining information on when and who is "on" is 
not straightforward even in a closed spatial domain. 

To resolve this difficulty, we first convert this local medium 
of a transmitting (or sensing) node from stochastic to deter- 
ministic using the equal sensing resolution concept in Section 
IIII-AI The medium here is thought of as a kind of two- 
dimensional space. The infinite plane is then regarded as the 
union, where this local medium of the same size is filled. 
Once the spatial boundary of this local medium is given 
as deterministic, N, the number of deployed nodes in this 
boundary is obtained using PPP. N a , the number of active 
nodes in this boundary should be found based on this N, and 
the transmission power distribution in this boundary can be 
found based on N and N a . 

These parameters are all analyzed based on the local shared 
medium called basic independent unit in this paper. By the 
homogeneity of PPP (Definition 2. in Appendix), the number 
of deployed nodes in the same area has the same probability 
mass function (PMF) regardless of its shape and location. 
Therefore, wherever the basic independent unit is, the number 
of nodes therein has identical PMF. Every related analysis is 
a function of this number of nodes. 



A. Carrier Sensing 

We introduce an effective CS range R so that a sensing node 
can sense any on-going transmission in this range and cannot 
sense any transmission outside of this area. By making the 
stochastic CS boolean, the complexity of the problem can be 
reduced. Then within the disk of radius -j, every node senses 
each other. We set this disk as a basic independent unit. 

CS is based on the threshold 7, i.e., if the sensed power level 
at a sensing node is greater or lower than 7, a sensing node 
regards the channel is busy or idle, respectively. Assuming 
there is only one interferer near the sensing node (i.e., we 
count only the nearest one if there are multiple interferers), 
the CS probability versus the distance to this interferer is 
calculated as followings: 
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Fig. 2. Pr[5i + v ^ 7] f° r varying distance d. pi ~ Exp(l/p). 7 is given 
as p/75 4 where p = O.OOl(iy) and noise v = 10~ 12 (W). In this example, 
R is calculated to be 70(m). 



P[Channel is busy] 



rPi 



(i) 



where pi is a RV representing the product of the fading effect 
and the constant transmission power from a typical node i, r 
is the distance between the sensing node and interferer, and 
v is the receiver noise power. Considering Rayleigh fading in 
d29i >. Pi follows Exp(l/p) with a constant transmission power 
P- 

Consequently, by the aforementioned definition of R, we 
convert the stochastic CS to a deterministic one. First, the 
average sensing area is calculated by integrating the parts of 
the circumference, of which the radius and the center are r 
and the sensing node, respectively. The CS probability of a 
point on this circumference is from (Q~|i. 



2irr ■ . 



2ixre 



{^7 + v> 7]dr 
r 4 



dr 



r 3/2 



t-N 



(2) 



Using the fact that the deterministic CS region should have 
the same average sensing area (sensing resolution) as the 
stochastic CS, the following equation is derived: 

7T 3 / 2 



-N 



= irR 



Finally, we get the CS distance R as follows: 



R 



1 / 7Tp 



y/2 \l-N 



1/4 



(3) 



(4) 



By this deterministic CS distance, the interference is regarded 
as boolean at a given distance rather than stochastic. 

Let us consider an infinite plane with the nodes randomly 
deployed. Delineating an arbitrary disk having of radius R/2 
in the plane, we denote this area by the sharing area^ where 

'This has the same object as the basic independent unit, i.e., the spatial 
unit where we can model the activity of nodes therein. However we use this 
term when explaining CS operation. 



every nodes in this area can sense other nodes' transmissions 
by the definition of CS range R. By PPR the number of 
deployed nodes in a sharing area follows Poisson distribution 
with the parameter Xtt(^) 2 as follows: 



[N = n] = 



{Mm 



.exp(-A^(|) 2 ) 



(5) 



Once we know N, we subsequently derive the number of 
active nodes and the power distribution at an arbitrary time 
instant in the sharing area. These will be explained in the 
following sections. 

B. Steady-state power distribution in a sharing area 



Definition 1. p on is the probability that there are on-going 
transmissions in a given sharing area at a certain time. 

So far we have derived F[N = n] and we will derive ¥[N a = 
a\N = n] based on p on . 

1) DCF dynamics in a sharing area: Consider the IEEE 
802.11 DCF protocol for CSMA/CA MAC. If all the trans- 
mitting nodes can sense each other in a sharing area and 
the given traffic to each node is a saturated one, we know 
the steady-state behavior in this area. As shown in [9] and 
subsequent research efforts, the backoff stage of each node 
in the network is random at a certain time and this can be 
elaborated through a two-dimensional Markov chain. We have 
two main quantities for addressing this: p c is the probability 
that the collision happens conditioned on the transmission of 
each node, and r is the transmission probability of a node at a 
randomly chosen time slot. These two quantities are derived by 
finding the steady-state solution of the discrete time Markov 
chain. 

By following the notations of [10], we have the BEB 
dynamics with m maximum backoff stage, K maximum retry 
limit and Wq initial window size. The backoff period durations 
are discrete uniform random variables 

given by: 



|W(0,2Wo-l), 
[U(0,2 m W - 1) 



0, . 

m 



,K-1 



(6) 



where U(a,b) is the uniform distribution with the outcome 
interval [a,b]. From this, the average backoff durations are: 



E[U {l 



(2W - l)/2 
(2 m Wo - l)/2 



for i = 0, . . . , 77 
for i = m + 1 , . 



,K-1 



(7) 

For a bigger size of data packet than an RTSthreshold in 
standard [11], K follows Long Retry Limit (the practical value 
is 4). Otherwise it follows Short Retry Limit (the practical 
value is 7^ Let 7r^ denote the relative frequency that a node 
enters the ith retransmission phase in the steady state. Then, 
it is shown that: 



7r, = (l-p c )p*(l-j£) 

-Refer to the simulation parameters in Table ITTT1 



(8) 
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and the probability r that a node transmits in a randomly 
chosen time slot is: 

(i - Pc)w Q (i - (2 Pc ) m ) 2 m w ( P ™~pK) ly 1 

2(l-pf)(l-2p c ) + 2(1 -pf) 2 
Again, p c is obtained through: 



i^TS 1 and GTS are the duration of RTS and CTS packet, 
respectively. 

T v has the PMF induced from (fTTT i such as p a (0), p a (l), 
and 1 — p a (0) — Pa(l) are f° r idle, successful transmission, 
and collision events, respectively. We derive the mean virtual 
time slot, E[T„] using this PMF in each mode. 

(9) ^BAS, 



Pc 



l-(l-r) 



(10) 



where N a is the number of active (contending) users. The 
probability that a transmitted packet encounters a collision, is 
the probability that, in a time slot, at least one of the remaining 
stations transmit. 

We can solve the system dynamics by solving two indepen- 
dent Equations (O and (TTOb and the existence of this solution 
is guaranteed by the fixed point theorem [9]. 

2) Power distribution in a sharing area: Since we know t, 
we can obtain the steady state power density. The probability 
that i nodes transmit simultaneously at an arbitrary time slot, 
given that N a transmitting nodes are deployed in a sharing 
area is: 

p a (m) = P[i = m\N a = a}= Mr m (l - r) a - m 

m = 0, . . . , a. (11) 

Each transmitting node's operation in a sharing area is 
synchronized since the medium is sensed perfectly and every 
node has the inner clock. Thus, idle time is also segmented 
into multiple slot times (a). Therefore, all events (idle time 
slot, successful and collision time slot) can be distinguished by 
their durations. At an arbitrary time, the sharing area medium 
is in one of three events and we call this random time slot 
as the virtual time slot. The virtual time slot has the random 
duration T v . We assume that the payload size is the same as 
PAY for all nodes. In basic mode, 



cr. 

for the idle slot time, 

T BAS {= PHY + [ (M ^ C + PAY) 1T S + sifs + ACK 
+difs), for the successful slot time, 
T C BAS (= PHY + \ {MAC + PAY ^ T S + difs), 
for the collision slot time, 



where PHY, sifs, ACK, difs are the durations for PHY 
header, SIFS time, ACK packet, DIFS time respectively. And 
M AC, R s and T s are the MAC header size, symbol rate and 
symbol duration respectively. Once the transmission starts, 
irrespective of success or not, the packet of size PHY + 
|- (mac+pay) i j,^ j s transm j tte( j fj rst A n( j m en tne remaining 

parts (sifs + ACK + difs or difs) are determined according 
to the existence of collision. Of course, in RTS-CTS mode, 
successful slot time and collision slot time will be changed 
into T S RTS (= RTS + CTS + PHY + [ (MAC + PAr) 1 T s + 3 * 
sifs + ACK + difs) and T RTS (= RTS + difs) respectively. 



= ap a (0)+T s BAS p a (l) 
E[T RTS ] 

= ap a (0)+T RTS p a (l) 



T C BAS (1- Pa (0)- Pa (l)), 
RTS (l-p a (0)-p a (l)). 



T 



The distribution of the number of concurrent transmis- 
sions (which is also the power distribution) is based on 
this PMF. In each virtual slot, the number of concurrent 
transmissions is various from to N a . And even when the 
number of concurrent transmitting nodes are given, the virtual 
slot consists of nobody transmission portion and somebody 
transmission portion. In the basic mode, nobody transmits 
during a. During sifs and difs in the both of successful 
and collision slot time, nobody transmits, too. During packet 
transmission time (PHY + \ ( MAC + PAY ) ]j; an d ACK) in 
a successful slot, one node transmits, while multiple nodes 
transmit during PHY + \ i MAC + PAY ) ~| m a collision 
slot. In RTS-CTS mode, the power density is changed in 
the same manner. Therefore, the actual power distribution, 
¥[j nodes transmit|iV a = a] is obtained using this distribution 
in the basic mode (Equation ( TTSi i) and RTS mode (Equation 
(fT3l)). The probability of busy channel in a sharing area p on 

■ v^q tjBAS/RTS, • \ -, ryBAS/RTS / a \ 
1S E i= l B a U) = 1 - B a (0). 



C. Number of active nodes in a sharing area 

Consider a given sharing area H. Let sensing area be 
defined as the sensing node's CS area excluding H. (See the 
asymmetric donut in Fig.[3]for the relation between the sensing 
and sharing areas.) The active node is defined here as the 
sensing node that has no on-going transmissions in its sensing 
area. In this subsection, we derive N a , the number of active 
nodes in H when the number of initially deployed sensing 
nodes in H is given The CS results for each sensing 
node in H is random. Therefore, N a is a RV that varies within 
[0,N]. 

For a sensing node in H to be active, two conditions must 
be satisfied: the CS result sensed from its sensing area must be 
idle and that sensed from H must also be idle. We investigate 
the former event in this section. The resulting N a becomes 
the input variable to the DCF behavior analysis in H, which 
is shown in Section IIII-B.ll 

In summary, since the areas of overlapped regions among 
sensing areas are varied due to the different locations of the 
sensing nodes in H and the activities of nodes in those over- 
lapped regions affect multiple corresponding sensing nodes 
simultaneously, the CS result of a typical sensing node in 
H cannot be modeled independently. Instead, by splitting 
a typical sensing area into three equally sized and disjoint 

3 By the homogeneous PPP property, it is not necessary to put any subscript 
or superscript of H on these variables. 
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(ap a (0) + (sifs + difs) Pa (l) + difs(l- Pa (0)-p a (l)), . ? =0 

B a AS U) = -p-L^T ' + r (MAC g +PAy) ir. + ACK) Pa (l), j = 1 (12) 

[v 1 {(PHY+\ ( MAC l P ^ ]T s )p a (J), 2<j<a. 

i (ap a {0) + [3sifs + difs)p a {l) + difs{l-p a (0)-p a (l)), j = 

B a TS U) = mT BT S] ■ \ ( RTS + CTS + PHY + \ iMAC + s PAY) ]T s + ACK)p a {\), j = 1 (13) 

[RTS-p a (j), 2<j<a. 




Sensing area of A 



Sensing area of B 



° Inactive node 
■ Active node 

X Center of 
sharing area 



Fig. 3. If we pick up a certain sharing area, every transmitting node in 
that area has its own sensing area. If there is no on-going transmission in its 
sensing area, that node is the active node. 



subspaces, we can quantize the area of the overlapped region 
in the unit of \H\ where | • | is the cardinality of the set. Then, 
by accounting for the combinations of the activities of each 
quantized area and the possible sensing node locations, we 
can find the distribution of N a as in (jT7j with the unknown 
variable p on . The following section contains the details of the 
derivations. 

1 ) Details of analysis: The following proposition gives the 
grounds for splitting a typical sensing area. 



Proposition 1. Each sensing area can be thought of as the 
union of three independent, equally sized and mutually disjoint 
areas, and the size of this area is the same as \H\. 
Proof. By the homogeneity of PPP, i.e., the point process is 
characterized by the area of the space not by the location or 
shape (See d24t in Appendix), the size of the sensing area is 
37vP , which is three times that of the sharing area. Therefore, 
if we measure the sensing area and compare it with the size 
of sharing area, we find that three sharing areas constitute one 
sensing area. ■ 

By Proposition 1, for one of the sensing nodes in H to be 
an active node, there should not be any transmission in any of 
these three independent areas, of which the size is \H\. Since, 
the point deployment is only specified by the size of the space, 
the probability that a node in H is an active node with only 
the sensing area considered is again a function of the activity 
in three subspaces, of which the size is the same as \H\. From 
this, that probability is 




Fig. 4. Sensing areas of multiple nodes of a sharing area. There are 
two nodes {A, B} in a sharing area. Sensing areas of two nodes include 
three black dots, {a, b, c} and {c, d, e} respectively. The hashed region is 
overlapped by both of sensing areas of node A and B. A dot (c) is included 
in this overlapped region, and means the quantized space occupied by this 
region. 



For example, if there is only one initially deployed node in 
H, the number of active nodes in that area has a simple PMF 
given by: 



P[JVa 



(1 -Pon) 3 , 
l-(l-Po 



a = 1 
a = 0. 



(15) 



(l-Pon) S 



(14) 



However, if multiple nodes are deployed in H, their CS 
operations are not independent of each sensing node since the 
sensing area of each node must partially overlap, as in Fig. |4] 

Due to this partial dependence of CS among multiple 
sensing nodes, the probability that a node in H is an active 
node is not straightforward like $15[ . The exact dependent 
CS probability cannot be derived since, to the best of our 
knowledge, even state-of-the-art research on stochastic geom- 
etry cannot reveal the exact distribution on the overlapped area 
when circles are deployed by PPP [12]. Therefore, we use 
a spatial quantization approximation, in which some evenly 
distributed dots are the measure of the overlapped sensing area. 
To do this, we introduce a large circle with radius that can 
cover the sensing areas of all possible nodes in H. We call 
this circle a big sensing area. There are eight dots in the large 
sensing area. They are evenly distributed in terms of the angle. 
As depicted in Fig. |5j any sensing node in H includes three 
dots, except for the case in which the sensing node is located 
exactly at the center of H. 



Proposition 2. When eight dots are evenly distributed on the 
circumference of a circle of radius R, the sensing area of any 
sensing node in H covers three consecutive dots. 
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Fig. 5. Each sensing area includes 3 consecutive dots in the big sensing 
area. Each dot represents one of three sensing areas. 




Fig. 6. The hashed region is the possible region of a sensing node having 
three-consecutive-dots {c, d, e}. This area has the angle of 4 which is split 
evenly by a direction to the center dot of {c, d, e}. If {c, d, e} are all 'off', 
this {c, d, e} becomes three-offs. The number of three-offs is bounded in [0, 8] 
in H. 



Proof. Let vector u characterize the shifting of the sensing 
area from the center of the large sensing area to a random 
position. The circumference of the original sensing area par- 
tially overlaps the shifted one. This partial circumference (the 
thick segment of the line in Fig. |5]l has an angle measured 
at the center of the sharing area, which is 2cos~ 1 (^-). Since 

< \u\ < % by the definition of u, the angle derived from u 
is bounded in the interval ,%]. This angle always includes 
three consecutive dots. ■ 

We measure the overlapped sensing area using these dots. The 
number of dots included in the intersection of any sensing 
areas denotes the corresponding overlapped area in the unit of 
\H\. Each dot represents the boolean indicating whether there 
is a transmission or not in the corresponding overlapped area 
in the unit of \H\. Each dot is "on" with p on and "off with 

1 — Pon, where p on for all dots is the same and the on-off 
boolean are independent of each other. The sensing area of 
any node in the sharing area can be mapped into one of three- 
consecutive-dots by Proposition 2. 

Conversely, the possible region of a sensing node having 
the same three-consecutive-dots exists in H as in Fig. [6] All 
sensing nodes located in this region follow the same CS result, 
i.e., on-off. For example, three consecutive dots of "off mean 
that the sensing area having these dots does not have any 
transmission in its area. The corresponding sensing node is 



Fig. 7. Example of counting three-offs. If the solid dot means 'off', r\ = 1 
in both cases, while D = 4 and 5 in (a) and (b) respectively. By counting all 
on-off combinations, we obtain Table U 
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28 
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256 



active. In contrast to this, if there is at least one of the 'on' dots 
among the three, this means that the corresponding sensing 
node is inactive. 

Any random position of a sensing node in H can be mapped 
into one of eight three-consecutive-dots and each dot has 
the probability to be on of p on . Therefore, the problem to 
find active nodes given N initially deployed nodes in H 
is solved in two stages. The first stage is to examine the 
case where three-consecutive-dots are all in "off states when 
each dot behaves on-off independently. The second stage 
is to find how many of these three-consecutive-dots in the 
"off state (simply say 'three-offs') are picked up through N 
independent choices. The number of all possible sequences 
of eight binary digits (representing the on-off states of a 
dot) are 2 8 = 256. By counting all combinations of these 
digits, we derive the PMF of the number of three-offs. Let 
us denote the number of three-offs by 77. For the values of 
7] = {0, 1,2,3, 4, 5, 6, 7, 8}, the number of occurrences of each 
r\ is {131, 56, 36, 16, 8, 8, 0, 0, 1}, respectively. These numbers 
are obtained simply by counting consecutive three-offs in all 
combinations of eight binary digits (or dots), as shown in Fig. 
□ 

These occurrences are also categorized according to the 
number of "on" dots among eight dots, which is denoted by 
D. By counting each combination, we obtain the quantity O®, 
the number of occurrences that D out of eight is "on" dots 
when there are ij of three-offs, as in Table U 

Therefore, the probability that there are r\ three-offs given 



x 



is: 



p n = F[i] three-offs] 

E°X(1-P™) 8 - D 



(16) 



D=0 



Once p on is given, p„ is determined using this equation. 
Since uniformly deployed sensing nodes in H have any three 
consecutive dots by Proposition 2, becoming an active node is 
equal to the event that the sensing area of that sensing node has 
one of 77 three-offs. Therefore, the probability that a sensing 
node uniformly deployed in H becomes an active node is | 
(rj = 0, 1, . . . , 8). A sensing node is located in H uniformly 
and independently. Hence, if there are n sensing nodes in H, 
the event that a out of n sensing nodes become active nodes 
is equal to the Bernoulli process with the parameter |. Then, 
N a ~ B(n, |). Let us define P n , a ,^ as the probability that a 
out of n sensing nodes have one of r\ three-offs when r\ is 
given. Then, 

Pn.a, n =P[N a = a\n,T]} 

n\ (rj\ a ( r]\ n ~ a 



Finally, N a is bounded in [0, N] and the conditional prob- 
ability that N a = a given N = n is: 



P[iV a = a\N = n] = J2 Pn,a,r,Pr,, 
r)=0 

a = 0, . . . , n. 



(17) 



D. Effective active node density 

We introduced a sharing area; using the spatial boundary 
of the boolean CS operation and derived the number of nodes 
distribution therein, which is P[N = n] as in (0, through the 
Section IIII-AI Using the sensing area concept, we derived the 
number of active nodes in a sharing area, P[N a ~ a\N = n] 
as in dPTl ) in Section IIII-CI Based on these derivations, the 
power distributions in a sharing area, (fT2l i and (13[ are derived 
in Sections IIII-B.ll and IIII-B.2I These results are all based on 
the value of p on , the probability that a sharing area is busy. 

The power distributions B a (j) of ( TTZb and ( fT~3b are con- 
ditioned on N a , while N a is conditioned on N, Therefore, 
the marginal PMF of p on is obtained by summing all of 
the probabilities on conditioned variables, as in ( TT~8b . The 
condition on N a is eliminated using ill) . The condition on N 
is eliminated by the homogeneity of PPP, i.e., every geometric 
subset satisfies Poisson distribution on the number of points 
in that subset, i.e., (0. 

In (II St . (a) is from the fact that the possibility of a 
sharing area being "on" is zero when there is no initial 
deployed nodes and active nodes, (b) is from (0 and ( TTTb . 
(c) is from ( fToT ). In this equation, p v includes p on term, 
which is also the marginal PMF from the unconditional 
event of channel busyness. Therefore, Equation ( TT~8T > has the 
unknown variable p on on both sides of the equation and 




Fig. 8. Each curve represents ELiE^n^.^Efl O^p^t 1 - 
Pon ) 8 ~ D } Ef=i Baij) of Equation [HQ when p on is varying and the value 
of n is given. The right hand side of Equation 081 is the linear combination 
of these curves. The intersection of this and y = x means the solution p on . 
This figure is for the case of A = 0.0004 and other parameters are in Table 

m 



the solution can be obtained by solving the weighted eighth 
order polynomial of p on , where the weight is the product 
form of Poisson and binomial probabilities. More specifically, 

Ea=l E„ Pn-aAED W„(l - Pon) 8 ^} £"=1 B a (j) is 

the eighth-order polynomial of the unknown p on and all other 



terms are known, and 



is a constant from 



Poisson distribution for the given n. Therefore, summing the 
product of these two terms for varying n from re = 1 to 
n = oo also results in an eighth-order polynomial. Since there 
is no general solution for polynomials with order higher than 
four and n goes to infinity, we cannot derive the closed form 
expression of p* n , the solution of p on . However, we can still 
get a solution by finding the intersection of the right hand 
side and left hand side of ( TT8I ). For computational simplicity, 
we use six instead of infinity for the value of n, as shown in 
Fig. [8] The rate of increase of the curves due to increasing n 
diminishes and the coefficient from Poisson distribution also 
rapidly decreases over higher n. Therefore the curve of n = 6 
is sufficiently close to the asymptotic curve. (We verify that 
this simplification is sufficient to reflect the original equation 
by comparing our analysis and simulation results in the next 
section.) 

If we obtain p* on , the distribution of the number of transmit- 
ting nodes in a sharing area can be derived as in (OH . where 
the number of transmitting nodes in a sharing area is denoted 
by Z. The expected number of transmitting nodes is derived 
from this result. 



\z 



(19) 



z=0 



The effective active node density is defined as the average 
number of transmitting nodes per unit area. Thus, we finally 
obtain the effective active node density as follows: 



A' = 



E[Z] 



(20) 



This is used in the cumulative distribution function (CDF) 
(f3Tb and the probability density function (PDF) (f32l > of the 
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(18) 



TABLE II 

Simulation Parameters for MHC and SSI processes. 



Fixed Parameters 



Transmission power (p) 


0.001(W) 


Background Plane 


Circle (B(0, R)) 


Radius of Background Plane (Rm 


= Rs) 


282 (m) 


Exclusion Ball Radius (r) 


70 (m) 


Number of Iterations (n) 


100000 


Channel Model 


where X ~ Exp(l/p) 


Varying 


Parameters 


Node Density (A) 




{1,2,3,4,5} * 10"* 



aggregate interference in the Appendix. We plot the resulting 
PDF for varying network parameters and compare this with 
the simulation results in Section [IV] 

IV. Verification of the analysis 

In this section, we first, plot the analyzed p on and A', and 
plot the PDF (T32t and CDF OTb of the aggregate interference 
using this A'. Next, we compare these derived results with 
the results of NS-2 simulations and two other dependent point 
process simulations. For all simulation scenarios, the exclusion 
radius r in dependent processes are given to 70 (m). In the 
NS-2 simulation, MAC/PHY parameters and channel model 
are given so that the CS radius R is determined to be between 
50 and 100 (m) for investigation of p on and A'. We later 
change these parameters for R to be equal to 70 (m) for the 
comparison with the dependent point process. 



A. Simulation setup 

1) Monte Carlo simulation for MHC and SSI: We deployed 
the points using MHC and SSI processes, explained in [5], 
with Matlab. The parameters used in the simulations are listed 
in Table [if]. For a given number of nodes, the aggregated 
power was measured at O the center of ball B with radius Rm 
or i?5 for MHC or SSI respectively. The number of deployed 
nodes was generated using CDF of Poisson distribution with 
the parameter X\B\ where A is the initial node density and 
\B\ is the area of ball B. We deployed these points uniformly 
and measured the aggregate power at O. We repeated this 
procedure for more than 100,000 iterations. 

4 Refer to the paper [5] for the meaning of parameters and details of the 
process. 



2) NS-2 simulation for PPP: To verify the analysis results, 
we conducted simulations using NS-2 [13]. Unlike the previ- 
ous version, NS-2 version 2.34 (released June 2009), includes 
wireless PHY and MAC layer patches for the realistic IEEE 
802.11 DCF standard [14]. This enabled us to realistically 
simulate the PHY and MAC stack of IEEE 802.11 DCF. The 
simulation parameters, which are the default ones for IEEE 
802.11a PHY and MAC, are listed in Table |ffl] They are from 
the previous research [15]. 

The overall simulation procedure consisted of genuine 
NS-2 simulations, pre-processing of the scenario, and post- 
processing of the data, as shown in Fig. [9] Generating 
the number of transmitting nodes using Poisson distribu- 
tion, deploying them uniformly, attaching designated receiv- 
ing nodes to each transmitting node, and generating traffic 
for each transmitting node were done in the pre-processing 
stage. For the saturated traffic situation, we obtained the 
time duration from maxijthe first transmission of node i} to 
minijthe last transmission of node i}. We call this the time 
window. Finding the time window, measuring the received 
power at the measuring node, recording the lasting time of 
each received power value, and accumulating all the measure- 
ments were done in the post-processing stage. 

In attaching the receivers to transmitters, we fixed the 
relative location of each receiver at 5 (m) right and 5 (m) 
up from its transmitter. To measure the aggregate power, we 
put the measuring node in the center of the simulation grid. 
This node then reported on the received power level and we 
recorded the lasting time and power level of each received 
signal. 

The simulation conducted in this paper is large-scaled; the 
reasons are two. First, each simulation per geometry scenario 
takes a long time. This time includes the simulation time in 
NS-2 and the post-processing time for handling received power 
instances. NS-2 traces all of the packet-level transactions 
with the received power recorded at every receiver. In post- 
processing stage, calculation of the received power from all 
of the on-going transmissions at a measuring node takes 
computation time. Moreover, the simulation time itself (not 
the computation time) has to be long enough to reflect the 
steady-state behavior, which theoretically requires infinite in- 
vestigation time. We consider at least 30 seconds per scenario 
as the simulation time (this takes about 2 h in real time using 
a Quad core i7 processor computer). Second, to get the sound 
statistical inference of PPP, we repeat per-geometry simulation 
many times. We do 50 repetitions, since all the resulting PDFs 
of aggregate interference obtained from the simulations show 
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TABLE III 
Simulation Parameters for NS-2. 



be saturated. 



Fixed Parameters 


Background Grid 


Regular Rectangular 


Grid Size 


500(m)*500(m) 


Transmission power 


0.001(W) 


Initial window size Wo 


16 


Maximum backoff stage m 


6 


CW min/max 


15/1023 


Slot Time (cr) 


9(us) 


SIFS 


16(us) 


DIFS 


sifs+2cr=34 (us) 


Short Retry Limit K 


7 


Long Retry Limit K 


4 


PLCP Preamble Duration 


16 (us) 


PLCP Header duration except Service field 


20(us) 


OFDM Symbol Duration 


4(us) 


IFQ Length 


50 


RTS MPDU + Service + Tail Field 


182(bits) 


CTS MPDU + Service + Tail Field 


134(bits) 


ACK MPDU + Service + Tail Field 


134(bits) 


Data Rate 


6(Mbps) 


Control Rate 


l(Mbps) 


Modulation 


BPSK 


Code Rate 


1/2 


Carrier Frequency 


5.18 (GHz) 


Preamble Capture Threshold 


2.5118 


Data Capture Threshold 


100 


Noise Floor 


10~ 12 (W) 


Data Type 


CBR 


CBR Rate 


5 * 10 B (packets per sec) 


Number of Packets in the application queue 


3000 


Varying Parameters 


Payload Size 


500 or 1000 (Bytes) 


RTS Threshold 


or 10000 


Node Density (A) 


{1,2,3,4,5} * 10~ 4 


Effective CS Range (R) 


50, 70, 100(m) 



the convergence before 30 repetitions. The simulation time for 
all 50 scenarios takes approximately four days on average. We 
repeat this process for each combination of PHY and MAC 
layer parameters. 

The saturated traffic was assigned to all transmitters so that 
there was no idle time by the traffic itself during the simu- 
lations. The background grid for all the simulation scenarios 
was always same: a square 500 (m) by 500 (m) in size. The 
transmission times for RTS, CTS, PPDU (PHY+MAC+PAY) 
with 500 B (or 1000 B) of payload and ACK were 52, 
44, 728 (or 1396), and 44 (ps), respectively (see Table [Hi] 
and [15]). We ignored the propagation delay, even though 
this exists in the simulator, since the value was quite small 
compared to the other transmission times. A generated traffic 
of 5,000,000 packets per second (pps) was given to all the 
transmitters. This means a value of 0.2 (/is) for inter-arrival 
time from the application layer to the physical layer, which 
is less than the whole transmission time for one successful 
packet, i.e., RTS + CTS + PHY + MAC + PAY + ACK + 
3 * sifs + difs = 948(/is), which is enough for the traffic to 



START 



Generate the random number N using Poisson 
distribution 

i 

Deploy H points (Tx) randomly in background 

plane 

1 

Attach points (Rx) to each Tx 

I 

Generate CBR traffic 

NS-7 simulation 

v 

Record the time window (TW): the first and the 
last time for saturated traffic 

4- 

Measure the received power and time at the 
measuring node during TW: PDF(i) 

•1 ~ 

Accumulate to the former one 

A(Q=A(f- 1) + /W(Q 

5 



Derive the final pdf,4 = 4^ 



END 



Fig. 9. Simulation flow. For a given combination of A, PHY and MAC 
parameter, this procedure is carried out. 



B. Discussions 




Fig. 10. P n versus A for all combinations of MAC layer parameters. 

1 ) Analysis results on p on and A': All p on in this subsection 
refer to p* on , the final solution of p on for simplicity of 
expression. In Fig. [TO] p on are shown for various combinations 
of MAC parameters. The factor that affects p on the most 
is the effective CS distance R, followed by A, the initial 
node density. As A increases, p on naturally increases due 
to the increased congestion level with semi-log style trend. 
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A' is higher than A" for varying A. Since A" is used for 
modeling dependent point processes, it cannot trace the real 
operation. As shown in the next section, our aggregate power 
distribution adopting this A' is the most elaborate among other 
point processes. Therefore, Fig. [TT] shows the gap in aggregate 
interference between simplified MHC and the real situation. 
This is significant because this simplified MHC expression is 
widely used in academia [3], [16], [17]. 

Finally, our new effective node density and the framework 
for deriving this can be applied to any other type of aggregate 
interference problem. In other words, whatever the appearance 
of the PDF, A, the original node density can be substituted 
for A', the new effective node density, using our derivation 
framework. 



Fig. 11. New node density A' for initial node density A and MAC layer 
parameter combinations. The thicker curves are using A" in {22) . 

If there is only reduced idle time factor without increasing 
collision probability, as in [8]'s ideal CSMA, the rate of 
increase of p on must be a log function style. However, our 
model reflects the increased collision probability and reduced 
idle time simultaneously, as obtained in real situations. The 
resulting p on shows a mixture of linear and log functions. 
Within the same R, the combination of mode and payload 
size that has the lowest p on is the RTS mode and short 
payload. In other words, the order of {RTS— 5005, Basic — 
5005, RTS - 1000B, Basic - 10005} is for lower p on . In 
general, RTS-CTS mode had a lower congestion level than 
the basic mode since the system cost was only the RTS-CTS 
packet collision and waiting time for retransmission compared 
to basic mode. A large payload in basic mode makes for 
a higher congestion level. However, in the cases of some A 
values, RTS-CTS packets were small enough to compensate 
for the payload size. Therefore Basic-500B might have higher 
p on than RTS-1000B. In this case, the effect of mode was 
weaker than that of the payload size. 

This p on was used in the new density A' d20l ) and we plotted 
this as shown in Fig. [TT] This figure shows that the smaller R 
makes for a higher A', which is the opposite case to p on - This 
is understandable since a smaller R signifies more insensitivity 
to the interference around. Therefore, we expected the result 
of R = 0, the ALOHA system, to approach the line A' = A 
in the figure. This line also represents wireless access systems 
that have no MAC, such as macro or femto cellular systems. 
The curve A' is a version of A filtered by the CSMA/CA and 
BEB mechanism. By showing the A' = A line and the curves 
together, Fig. [TT] also addresses how large the gap between 
these two node densities is and how effectively CSMA/CA 
MAC operates. The bold curves are from: 

A" = l-c*P(-^) (22) 

where D is the exclusion distance with an ambiguity. However, 
we put R into D in this figure. This expression is the 
approximated node density for modeling MHC adopted in 
[3] and [16]. As can be seen in the figure, the variation of 
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Fig. 12. Probability density of aggregate interference in the condition of 
RTS mode, 500B payload and R = 70(m). Other parameters are in Table III. 




Fig. 13. Cumulative distribution of aggregate interference in the condition 
of RTS mode, 500B payload and R = 70(m). 

2) Comparison of the resulting aggregate interference with 
the simulation: Using this A' and shot noise analysis, we 
plotted the distribution of aggregate interference. The PDF 
and CDF of the analysis in each node density showed high 
correlations with those of the NS-2 simulations, as seen in Fig. 
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Fig. 14. Probability density of aggregate interference when A is 0.0001. 



Fig. 16. Probability density of aggregate interference when A is 0.0005. 
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Fig. 15. Probability density of aggregate interference when A is 0.0003. 
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Fig. 17. Probability density of aggregate interference in the case of 500B 
payload and R = 70. 



n~2l and Qj] We depict all the PDFs in a log scale. Although 
at first glance they resemble a log-normal distribution, they 
are asymmetric based on the main lobe. They show heavy 
tail properties even in a log scale axis by investigating the 
following equations: 
Normal scale: 



lim e at Pi[I > t] 
Log scale: 



lim e at Pr[I > logt] = lim e° 



lim e at 1 - erfc 



1 — erfc 



4y^ 



where / is RV of aggregate interference. 

Therefore, they are definitely neither normal nor log-normal 
distributions. This is notable as some research efforts in the 
signal processing field assume that the aggregate interference 
follows normal (in dBm unit) or log-normal (in W unit) 
distributions. For the other features, the higher the mean of 
the aggregate power, the lower the probability of that mean 
value. Therefore, low-mean high-probability and high-mean 
low-probability patterns are shown in all the results. This is 
because the total sum of the probability is fixed to 1 and the 
x-axis is log-scaled and not linear. 

Compared with dependent point processes, at any given 
A value, our analysis is the closest one to the simulation 
results, as depicted in Fig. [l4lFig. [161 The interference of 
MHC is always less than that of SSI since MHC is the lower 



bound of SSI, which is commented on in [5]. However, these 
two are not as sensitive as our model to variations in node 
density. Moreover, these two do not have sufficient MAC and 
PHY layer parameters to reflect the real situation, while our 
analysis can model any combination of the system parameters, 
as in Fig. [T7] As shown in the figures so far, our model 
of interference has slightly lower values than that of the 
simulation in each case. This is mainly because the simulator 
allows the capture situation. In our analysis, the collision 
between transmitters is regarded as a failure of transmission 
and this increases each collided nodes' backoff stages. In 
contrast, there might be a successful transmission even when 
multiple nodes in a CS area are transmitting at the same time. 
This is because if the ratio of one incoming signal to the 
others is guaranteed to be higher than a certain threshold, 
this stronger incoming signal can always be decodecd. In our 
analysis, we ignored that situation in order to simplify the 
analysis. Nevertheless, our result is competitive. 

From these results, we learned the following lessons: If the 
network is required to maintain a lower interference than a 
certain level, there are multiple combinations of parameters 
that need to be controlled. Since the PDF of the aggregate 
interference is a function of A', there are multiple combinations 
of parameters that can result in the same value of A'. Those 

5 The patch in NS-2 has two capture thresholds which are the preamble 
capture and the data capture. Thus, these two capture events can also occur 
in the simulations. 
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controllable parameters are R, transmission mode, payload 
size, etc. This can be used for the interference management 
in uncontrolled interference limited systems such as cognitive 
radio networks. In addition, this statistical information helps to 
quantify the spectrum opportunity in cognitive radio networks 
such as exploited in [18]. In addition, the interaction among 
heterogeneous systems using the same frequency band such as 
in WLANs, ZigBee, and Bluetooth can be understood using 
this analysis from the perspective of the mutual interference. 
This knowledge can also be used to measure the energy level 
in the air when wireless power transfer (WPT) systems have 
some conditional activation mechanisms. 

V. Conclusion 

In this paper, we analyzed the aggregate interference from 
randomly deployed CSMA/CA nodes. CSMA/CA MAC op- 
eration on the stochastic geometry makes the understanding 
of a node's behavior difficult. However, the homogeneity of 
PPP provides a clue to solving the problem. We derived 
the effective active node density by spatially quantizing the 
infinite space and analyzing the steady-state power distribution 
on this quantization unit. Verification is also nontrivial since 
repetition is needed in order to get sound statistical inference 
from random topology. Although the elaborate closed form 
expression of the interference distribution cannot be obtained, 
the analysis framework proposed in this paper shows high 
applicability to many related problems in modern wireless 
networks where the amount of interference in a random node 
geometry significantly affects system performance. 

Appendix 

This section can be replaced with reference [6] as long as 
the readability does not deteriorate. 

A. Aggregate power from Nodes under the PPP Geometry 

We consider the infinite planar where the transmitting nodes 
are deployed randomly at positions specified by a Poisson 
distribution. For the general dimension space M. d and the 
locally finite non-null measure A, the PPP can be defined: 



Definition 2. The PPP $ of intensity measure A is defined by 
means of its finite-dimensional distributions 



P{$(Ai) = m, . . . , $(Ak) = n k } = ]J 



i=i 



(23) 

for every k — 1, 2, . . . and all bounded, mutually disjoint sets 
Ai for i = 1, . . . , k. If A(dx) = Xdx is a multiple of Lebesgue 
measure (volume) in M d , we call $ a homogenous PPP and A 
is its intensity parameter. 

Simply, in 2D homogeneous PPP, the probability that the 
area S has k nodes is: 



P(k) = e 



k\ ' 



Proposition 3. (Proposition 1.2.2 in [6]) The Laplace func 
tional of PPP of intensity measure A is 

S^ d (l-e-f^)h(dx) 



(25) 



(24) 



where A is constant for all disjointed and bounded area S. 



where d is the dimension of point deployment and / is the 
probability function of the point position. See the [6] for the 
proof. 

A shot-noise field is non-negative vector random field 
defined for all y in some Euclidean space and which is a 
functional of a marked point process <P. Here is a list of 
arguments dimensions related to this field: All y £ M. d ; 
I$(y) £ K +fc for all y; and this is generated by a marked 
point process $ = ^ st Xi>m A on R d with marks in M. 1 where 
£ is a Dirac measure; and non-negative response function 
L = (L±, ...,Lk). Therefore the total dimension input and 
output arguments of shot noise field is R d xR d xR l — s- 



Definition 3. The shot noise field associated with the 
marked point process <3? and the response function L = 
(Li,L 2 , ...,L k ) is 

hiv) = ( I i(y),---,Ik{y)) = ^2 L{y,x u mi), (26) 

where the sum is evaluated component-wise for the vector 
response model L. 

Proposition 4.(Proposition 2.2.4 in [6]) Suppose that $ is an 
independent marked PPP with intensity measure A and mark 
distribution F x (dm). Consider the Shot Noise I(y) — I$(y) 
with response function L. Then 

= exp(-/ / (l - e -^=i SiLi to' x < m A F x (dm)A(dm)\ .(27) 

See the [6] for proof. 

Here Li(y,x,m), the element of the response function L 
represents the i-th channel link gain function out of k ones 
with the measuring location y, transmitter position x and 
mark to. The total power received from this collection of 
transmitters at a given measuring location is in essence a 
shot-noise field. Applying this proposition to our CSMA/CA 
random network in an infinite planar with omni-directional 
antennas, the shot noise model consists of 4 components. 1) 
at a measuring location in 2-D planar (df = 2); 2) collection of 
points representing the locations of transmitters on the plane 
R 2 (d — 2); 3) marks m, = pi representing the powers of 
the transmitters (£ = 1); and 4) a scalar (k = 1) response 
function L(y,Xi,pi) = Pi/l{\xi — y\), where l(r) is the 
omni-directional path-loss function of distance r. According to 
whether the points deployment and the assignment of mark to 
a point are deterministic or stochastic, the standard stochastic 
scenario is categorized. Our model is one of the four 

scenarios, which means the mixture of Poisson point process 
and exponential RV assignment for a mark assignment. 

6 M/M is the expression of Kendall like notation which means Markov point 
deployment and Markov mark assignment in the standard stochastic scenario. 



14 



For the general l(r) and G(s) — Pi[pi < s], Proposition 4 
leads 



£/(„)(«) = Ci(s) 



cxp|-2^A^ r(l - C p (s/l(r)))dr^ , 



(28) 



where C p (t) = / n °° e~ ts G(ds) is the Laplace transform of the 
transmitted power. 

We use the simplified pathloss function, l(r) = (Ar) 13 
where A > 0, (3 > 2 for an omni-directional antenna case. For 
the power reception model, we use the virtual power model 
where each transmitting node assumes to emit a random power 
not the deterministic one and there is only pathloss effect in 
the channel link gain. In other words, the transmission 
power at node i equals PiF where Pi is the constant value 
and F is independent and identically distributed RV. And the 
received power at the measuring location originated from the 
node i is l(ri)pi where r,; is the distance between node i 
and the measuring location. For modeling Rayleigh fading 
channel, we assume that F is exponential RV with unit mean. 
And Pi is constant p for all transmitting nodes i. Then, pi 
is exponential RV with the mean p ■ 1, i.e., pi ~ Exp(i) and 
G(s) = Pr[pi < s] = l — e~p and we get following equations, 
which are Eq. (2.25) and (2.26) in [6]: 



Ci(s) = cxp 



(29) 



27rr(2//3)r(l-2//3) 



where K(f3) = p J^ /p) 

Unfortunately, there is no known closed form of I for 
general j3 but 4 since the inverse Laplace transform does 
not exist for all j3 cases. From now on, we only consider 
A = 1 and /3 = 4 case. Let / and F denote PDF and CDF 



_ dF(t) 
dt 



holds. From 



of the same RV respectively, then f(t) 

Cf(s) = s ■ Cf{s), we have the Laplace transform of CDF of 
RV I, _ 
cxp{-AX(4) v ^} 
C Fl (s) = v - . (30) 



[3] F. baccelli, B. Blaszczyszyn, and P. Muhlethaler, "An Aloha protocol 

for multihop mobile wireless networks, " IEEE Trans, on Information 

Theory, Vol. 52, No. 2, pp. 421-436, 2006. 
[4] J. Andrews, R. K. Ganti, M. Haenggi, N. Jindal, and S. Weber, "A 

primer on spatial modeling and analysis in wireless networks," IEEE 

comm. magazine, Vol. 48, No. 11, pp. 156-163, 2010. 
[5] A. Busson, G. Chelius, and J.-M. Gorce, "Interference Modeling in 

CSMA Multi-Hop Wireless Networks," Technical report INRIA, inria- 

00316029-ver3, pp. 1-21, 2009. 
[6] F. baccelli, and B. Blaszczyszyn, "Stochastic Geometry and Wireless 

Networks Volume 1 : THEORY," Foundations and Trends in Networking, 

Vol. 3, No. 3-4, pp. 249-449. 2009. 
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Using the fact that inverse Laplace transform of ex P( j s 
erfc(^) when a > 0, we obtain CDF and PDF of RV I 
similarly in [10]: 



Fi(t) 



//(*) 



erfc 



exp( : 



Att 2 ^ 

4Vt 

-^ 4 A 2 p V 
16* /' 



r 3/2A 

P 



4(1)3/2 



(31) 



(32) 



where erfc(a;) 
function. 



2 



is a complementary error 
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